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I. INTRODUCTION 

Though there is no evidence for the presence of a fundamental scalar field, it has played 
an important role for the study of the earliest stages of the universe |l|] ||^. According to 
this interest, the boson star, which was first discovered theoretically by Kaup and by 
Rufiini and Bonazzola [Q, has become an interesting subject in its own right as well as in 
the problem of the cosmological missing mass problem. This configuration has been studied 
in various theories [0 B- Their stability |T^ and structure [|rT| have been also 



investigated in some literatures, (for reviewing the boson star, see Refs. |T^ and [13], and 
other references are cited therein). 

Boson stars are formed by soliton-type configurations which are held together simply by 
their self-generated gravitational field and are only prevented from gravitational collapse by 
the Heisenberg uncertainty principle. Thus, we can easily calculate the order of radius Rq, 
maximum mass M^"-^ , and energy density po of a spherically symmetric stable boson star 
lias 

i?o~l/m, Mo^'^^-MpVm, po ~ M^m^ (1.1) 

where m and Mp are the boson mass and the Planck mass, respectively (we use the unit set 
hy h = c = 1). Then, the magnitude of the central value of the scalar field leading to the 
most massive stable boson stars has the order of the Planck mass, |0| ~ Mp. In Eq. (|l.l| ), 
we see that the maximum mass of a boson star is much less than the Chandrasekhar mass, 
Mch ~ Mp/m?. For a bosonic particle with mass 100 GeV, comparing with a neutron 
star, the maximum mass, radius, and density of a boson star are 10^% ~ 10"'^^ Mneutr on, 
lO^^^m ~ 10'^"^ Rneutron, and lO^^kg/m^ ~ lO^Vneuir-on, respectively |13| . 



The maximum mass of such a stable "mini-boson star" is too small to be a candidate for 
the missing mass. It is worthy of notice that why a boson star has a much smaller mass than a 
typical fermionic star is due to the fact that a boson star is only prevented from gravitational 
collapse by the Heisenberg uncertainty principle instead of the Pauli exclusion principle. In 



other words, the Heisenberg uncertainty principle is characterized by a much smaller length 
scale than that of the Pauli exclusion principle. It has been enhanced by Colpi, Shapiro 
and Wasserman ^, who considered the boson star formed with self- interacting scalar fields. 
They have shown that when the order of the coupling constant A is about unity, the boson 
star may have the maximum mass comparable to that of a fermionic star. 

^max _ o.22A^/2 MpVm ~ 0.22X^^^M^h, (1-2) 

where A is a dimensionless quantity defined by A = XAIp/Airm'^. The origin of this drastic 
change is that the introduced self-interaction term in matter Lagrangian generates a repulsive 
force. Thus, the characteristic length scale of the theory becomes large, and the boson star 
becomes large and massive. 

In this paper, we shall make analytic derivation for the maximum mass in Eq.( [L.2D 
calculated by numerical analysis. (Analytic derivation of maximum mass of boson star was 
also appeared in Ref. fl^.) Our analytic derivation will be also applied to the boson star 
formed by scalar fields with any higher order self-interactions. Then, a general form of the 
maximum mass including the self-interaction effect is obtained and it will be shown that 
the contribution of the higher order self-interaction terms to the maximum mass decreases 
as {m/MpY power. In order to confirm our analytic calculations, we shall make numerical 
analysis for the scalar field with the quartic and the sixth order self-interactionsQ. In Ref. 
|T3| , the authors chose the magnitude of the central value of the scalar field to be the order 



of the Planck mass, |0| ~ Mp, even for the self-interacting scalar fields. However, we shall 
show that when the contribution of the self- interaction becomes important, it has the order 
of the mass of the scalar field, |0| ~ mA~^/^ (This relation seems to be singular in the limit 
of A -^ 0. But, our derivation of the relation will be given only in the case that the limit is 
invalid.) This is consistent with the numerical calculation given in Ref. |^. 



^For the boson star formed by the scalar field with the sixth order self-interaction, see Refs. |15| 
and lie I 



This paper is organized as follows; In Sect. II, the order of the maximum mass of a stable 
boson star formed by the scalar field with the quartic and the sixth order self-interactions is 
analytically derived. Then, the formulation is generalized to the scalar field with any higher 
order self-interactions. Our derivation is confirmed in numerical analysis given in Sect. III. 
In Sect. IV, we briefly summarize and discuss our results. 

II. ANALYTIC DERIVATION FOR MAXIMUM MASSES 

In this section we make analytic derivations of the maximum masses of boson stars 
formed by self-interacting scalar fields. Firstly, this is done for the cases that the matter 
Lagrangian includes the quartic and the sixth order self-interaction terms. Then, we are 
going to extend our formulation for any higher order self-interacting bosons. 

We begin with the Lagrangian of the scalar field minimally coupled to gravity given 
by 

^ = -^^^>;0;. - ^^'101' - \m" - ^7101', (2.1) 

where A and 7 are the coupling constants. 

As mentioned above, when one is not concerned about the self-interaction, the magnitude 
of the central value of the scalar field leading to the most massive stable boson stars |0| is the 
order of the Planck mass. When the contribution of the self- interaction becomes important, 
putting 101 is of the order of the Planck mass, however, will give rise to some inconsistencies. 
This is because the self-interaction generates a repulsive force and the characteristic length 
scale of the theory becomes larger. To see this, consider only the quartic self-interaction 
term (7 = 0) in the Lagrangian (|2.1|) . If one set the order of |0| with that of the Planck 
mass, the energy density including the quartic self-interaction term p\ (Eq.(14) in Ref. |I3[ ) 
becomes 

Pa ~ m2|0|2 + A|0|' ~ m'Mlil + A). (2.2) 



Comparing Eqs.( p..l|) and ( |2.2| ), we see that this corresponds to the energy density of a star 
formed from non- interacting bosons with rescaled mass m -^ m(l + A)^/^ and the order of 
the radius of the boson star is given by 

Rx--{l + K)-^'Vm. (2.3) 

Thus, the maximum mass of the boson star becomes 

M];^"^ ~ PxRl ~ (1 + ky^l^Mllm ~ K~^'^Ml/ni. (2.4) 

Note that the maximum mass in (p.4|) is not equivalent to the numerical result (|1 . 2|) given 
in Ref. 0. Moreover, p\ in Eq.(|2.2|) is much greater than po in Eq. (pTTP and R\ in Eq.( p.3|) 
much less than Rq in Eq. (|l.l| ). It seems to be inconsistent with the fact that the introduced 
self-interacting potential generates a repulsive force. 

Now, turn to the energy density including the quartic self-interaction term. Since we 
expect from the field equation that the contribution of the self-interaction term is comparable 
to the mass term in the energy density, m^|0p ~ A|0|^, the order of \(f)\ becomes 



ni 



2UI2 



^T7?~^(l)' 101 ~ MpA-i/2 ^ mA-i/2_ (2.5) 



"M2 



V 

Note that \(f)\ has the order of the scalar mass not the Plank mass. In fact, we can see that 
in Fig.l of Ref. |^, the central value of |0| decreases with increasing A. (The Fig.l of Ref. 
^ was plotted only up to the value A = 300. However, the dimensionless constant A is very 
large, A k, 10'^^ ~ 10'^^. For this value, the central density ctq in the figure would be very 
small.) On the other hand, the Eq.( p.5| ) seems to be singular in the limit of A ^ 0. However, 
as mentioned above, the Eq.( p.5|) is only available in the case that the self-interaction term 
is so large to be comparable to the mass term in the energy density. The energy density 
becomes 

PA~"^'|0|' + A|0|'~m2M2A-\ (2.6) 

Again, this corresponds to the energy density of a star formed from non-interacting bosons 
with rescaled mass m — > mAr^f'^. Thus, the radius of the boson star becomes 



Rx ~ A'/ym. (2.7) 

In Eqs.(p.6D and (p.7|), the density and the radius become dilute and large, respectively, and 
it is consistent with the fact that a repulsive force is involved. Finally, using Eqs.( [^.6D and 
( [^.7D , we obtain the maximum mass of the boson star equivalent to the numerical result in 
( |1.2D up to constant 

j^rnax _ f^y^M^/m ~ X^^^M^h- (2.8) 

Let us now consider the sixth order self-interacting bosons. Following the above formu- 
lation, the central value of the scalar field is given by 

A|0|^ + 7|0r f\<P\V 



m2|0|2 \Mp^ 



Af^Vn-r 



4 



1^1 (A2 + r)~C(l), (2.9) 



Mp 



\4>\^M,{A' + Tr'/\ (2.10) 

where F = 7Mp/(47rm)^. Note that we used the mean field approximation on the second 

step in Eq.(p.9D, i.e., the relation |0| ~ MpA~^^'^ given in Eq.p.5|). Thus, the energy density 
and radius of the boson star are given by 

p,, ~ m'\4>f + A|</)|^ + 7|</>r ~ m'M'^iA' + T)-'/', (2.11) 

Rx^r^iA^ + ry/ym, (2.12) 
respectively. As a result, the maximum mass becomes 



M— ~ (A^ + ry/^M^^/m r.(x'+^ (^) 1 



V4 / / X 2\ V4 



m 



K/^ -\>^ +'^[-mJ I ^^'^' 



(2.13) 



where 7 = F(47rm)7Mp2 = 'jMj,. If 7 = 0, the Eq.( p33D becomes Eq.([T3). Note that from 



the field theoretic point of view, 7 is about unity and the contribution of the sixth order 
self-interaction to the maximum mass is minor as the order of the inverse square Plank mass 
comparing with the quartic self-interaction term. 



Our above analytic derivations can be easily extended to the case that the matter La- 
grangian involves any higher order self-interaction terms. The self-interaction terms can be 
expressed in terms of dimensionless parameters X{2k+2) as[| 



A 



(4)1 



|4 + MU|6 + MU|8 + . . . + ^(2«+2) ,^,2n+2 



M2 
p 



p 



M2n~2 
p 



(2.14) 



These higher order nonrenormalizable terms are naturally expected from the gravity as low 
energy effective theory and the dimensionful constants are given as the power of \{2k+2) / M^ . 
Then, the magnitude of the central value of the scalar field becomes 



|(2n+2) ~ m 



Sfc=iA(2A:+2) 



\ 2(n-fc) / \ 2fc-2' 

m \ I m \ 



l(2n), 



M„ 



-l/2ra 



(2.15) 



where |0|(a;) is the magnitude of the scalar field considering the order of the self- interaction 
to be k. The radius and the maximum mass are given by 



R 



(2n+2) ~ 2 

m 



SLiA 



\ 2{n-k) / \ 2fc-2' 

m \ m \ 



(2fc+2) 



l(2n). 



M„ 



l/2n 



(2.16) 



M, 



(2n+2) 



Sfc=iA(2A:+2) 



\ 2(n-k) / \ 2fc-2' 

m \ m \ 



l(2n). 



M„ 



l/2n 



M,h, 



respectively. The Eq.( |2.17| ) can be rewritten by 



M, 



(2n+2) 






+ • • • +A 



(2n+2) 



m 

m: 



(2n-2) 



+ cross terms 



l/2n 



M, 



ch- 



(2.17) 



(2.18) 



As a result, in Eq. (|2.18|) , the contribution of the higher order self-interaction terms to the 
maximum mass decreases as {m/MpY power. 



The higher order self-interaction potentials arise from effective theories is discussed in Ref. |16|. 



III. NUMERICAL ANALYSIS FOR MAXIMUM MASSES 

In this section, we make numerical analysis to conform maximum masses of stable boson 
stars calculated in previous section. We only consider up to the sixth order self-interaction 
term. The equations of motion generated from the Lagrangian ( pj.l| ) are given by 



G^^ = MGTj^, 



(3.1) 



V^V^ct) - m^ct) - A|0| V - 7I0IV = 0, 



(3.2) 



where the energy-momentum tensor T^y is given by 



Tji = Ig''' {<p:a,. + </';<x</'r.) - l^i^ 



'g^-<P*<j,.^^ + m'\<pf + lx\<P\' + y<pf 



(3.3) 



2 '"' 3 

We seek a spherically symmetric, time-independent solution represented by the line element 
in Schwarzschild coordinates 



ds' = -B ir) dt' + A (r) dr' + r'dVt. 



(3.4) 



Then, we may choose the ansatz for the scalar field </> as 



(r, t) = $ (r) e 



-iujt 



(3.5) 



and the Eqs. (|3.1| ) and (|3.2|) are written by 



A' 1 / 1 



fi+^r +2" +3^ +^4-' 



(3.6) 



x^\ A 



ABx 



5/23 A 



/\2 



(3.7) 



// /^2 B' A'\ , ^ 

a" + - H -\a' + A 

\x 2B 2A 



^-l\a-Aa^- ra^ 



where prime denotes d/dx, and x = mr, a = (AnG)'-^ $, and Q = uo/ra. If we write 



Aix) 



1- 



2M{x) 



X 



(3.8) 



(3.9) 



the Eq. (|3.6|) can be rewritten by 



M' = x^ 



2\B I 4 6 2A 

The total mass of a boson star is given by 



(3.10) 



Mx^ = M (oo) (Ml/m) . (3.11) 

We require that boundary conditions be A (0) = I, cr (0) = (Jq, o' (0) = and o (oo) = so 
that our solution becomes a nonsingular and asymptotically flat one with finite mass. 

First of all, we consider A = case. We can see in Fig.|l] that the boson star mass 
increases with increasing F. Explicitly, from Fig.^] plotted the maximum mass of boson star 
as a function of F, we find out that the maximum mass of boson star behaves as 

/ \ 1/2 

j^vaax ^ o.24Fi/^MpVm ^ 0.247^/^ ^ M^h- (3.12) 



This result corresponds to the maximum mass in ( p.l3| ) with A = 



Let us now consider A 7^ case. Numerical results of A4 (00) for A = 10, 30, 100 are 
shown in Figs.^,^^. We can see that F dependence is not sensible for large value of A. From 
Fig.^ in which we plot the maximum mass of boson star as a function of F for several different 
A values. It can be seen that dotted lines representing numerical results approach the solid 
lines drawing (0.22^A^ + 0.24^F) . Thus, the maximum mass formed by the self-interacting 
scalar field is given by 

j^rnax ^ |^o.22^A2 + 0.24^f) ^^^ M^/m 

0.22^ A^ + 0.24^7 ("^M M,,. (3.13) 

This is equivalent to Eq.( |2.13| ) up to constant factors. 

Figs. 1,3,4, 5 show that the numerical noises become noticeable when F reaches about 
10000. However, the maximum masses of boson stars are not affected by this part of the 
numerical error. 



IV. SUMMARY AND DISCUSSIONS 

Boson stars are only prevented from gravitational collapse by the Heisenberg uncertainty 
principle whose characteristic length scale is enormously small than that of the Pauli ex- 
clusion principle. Thus, the radius and mass of a boson star are much less than those of a 
fermionic star. However, Colpi, Shapiro, and Wasserman have shown that introducing the 
self-interacting scalar fields, the maximum mass of a stable boson star may be comparable 
to the Chandrasekhar mass. 

In this paper, we have made analytic derivation for the maximum mass and have shown 
that it is equivalent to their numerical result. In addition, it has been shown that the 
radius and density of the boson star formed by the self-interacting scalar fields become 
larger and dilute compared with those of the boson star formed by the scalar fields without 
self-interaction. The origin of this drastic change is that the introduced self-interaction term 
in matter Lagrangian generates a repulsive force. Thus, the characteristic length scale of the 
theory becomes large, and the boson star becomes large and massive. We also generalized 
our analytic derivation for the boson star formed by scalar fields with any higher order self- 
interactions. It has been shown that the contribution of the higher order self-interaction 
terms to the maximum mass decreases as {m/MpY power as expected in the viewpoint of 
the field theory. Our analytic calculations have been confirmed in numerical analysis in 
which we considered the scalar field with the quartic and the sixth order self-interactions. 

Recently, rotating boson stars have been also considered [|l7l. This generalization is likely 
to be natural and it would be interesting to extend our analytic calculation to the rotating 
boson stars. 
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FIGURES 
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FIG. 1. Boson star mass M^ as a function of ctq when A = 
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FIG. 2. Maximum mass of boson star as a function of F. 
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FIG. 3. Boson star mass Mx^ as a function of o"o when A = 10 
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FIG. 4. Boson star mass Mx^ as a function of uq when A = 30 
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FIG. 5. Boson star mass M\^ as a function of ctq when A = 100 
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FIG. 6. Maximum mass of boson star as a function of F under A = 0, 100, 200, 300. 
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